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Abstract. We study theoretically the enhancement of spontaneous emission in wire 
metamaterials. We analyze the dependence of the Purcell factor dependence on wire 
dielectric constant for both electric and magnetic dipole sources, and find an optimal 
value of the dielectric constant for maximizing the Purcell factor for the electric dipole. 
We obtain analytical expressions for the Purcell factor and also provide estimates 
for the Purcell factor in realistic structures operating in both microwave and optical 
spectral range. 
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1. Introduction 

Wire metamaterials are composed of arrays of optically thin metallic rods embedded in a 
dielectric matrix [1] . Experimental realizations of such structures span from microwaves 
[2] to optics [HI m ini El [7], and they are very promising for a number of applications 
including the subwalength transmission of images [2] , negative refraction phenomena [8] , 
superlensing [9], and biosensing applications [5]. Strong enhancement of the Vavilov- 
Cherenkov radiation and peculiar dipole emission patterns have been also predicted for 
wire metamaterials [TT| [T2] . 

Specific property of wire metamaterials is the strong spatial dispersion of the 
effective dielectric constant [13], manifested as the formation of so-called TEM modes in 
addition to the TE and TM modes of an ordinary uniaxial structure. Due to a finite value 
of the dielectric constant of wires, TEM modes acquire hyperbolic dispersion [Tl]. In this 
regard the wire medium represents a particular class of metamaterials with hyperbolic 
isofrequency surfaces, which have recently attracted a lot of attention [I5l HE]. The 
specific feature of the hyperbolic metamaterials is diverging photonic density of states, 
promoting high spontaneous decay rate of the embedded light source [T7l [18]. In wire 
metamaterials this effect should lead to large Purcell factor as well. 

Although the Purcell factor in hyperbolic metamaterials stays finite due to the 
discreteness of the actual structure [19], [201 1211 122] , its value is strongly sensitive to the 
structure geometry. Consequently, the detailed theory recently developed for planar 
metal-dielectric metamaterials [13 [IHl [20] is not applicable directly to nanowire arrays 
and novel study is required. The crude estimation of the Purcell factor is provided by the 
density of states enhancement of TEM modes~ (A/a)^ [23], where A is the wavelength 
and a is the structure period. Our goal is to perform more comprehensive analysis which 
takes into account finite dielectric constant of the wires and arbitrary spatial position 
of the source in the metamaterial unit cell. We consider both electric and magnetic 
dipole sources, because high sensitivity of the enhancement factor to the dipole source 
is known already since the original work of Purcell [24] [25] . 

The paper is organized as follows. Section [2] presents our theoretical model. 
Section [3] is devoted to the analysis of the eigenmode dispersion of the structure. 
Section [H contains numerical and analytical results for the Purcell factor. Discussion 
of the Purcell factor attainable in different experimental conditions is summarized in 
Sec. [5l whereas the last section concludes the paper. 

2. Local field approach 

The structure under consideration is illustrated schematically in Fig. [H It consists of 
identical infinitely long wires with radius i?, arranged in infinite square lattice with 
period a and parallel to z axis. Wires are characterized with dielectric constant ^wire- 
Electric (p) or magnetic (m) dipole source is embedded in the structure at the point 
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ro = pQ + zqz. Electric field satisfies the following equation: 
V X V X - q^e{p)E = Airq^P 

where 

e{p) = 1 + (£„ire - 1) ^ e{R -\p- p^\) 

j 

is the dielectric function (here 9 is the Heaviside step function), 

p6{r - ro), 



(1) 
(2) 



Q 



-V X m6{r — ro). 



(3) 



is the polarization term, describing the point electric or magnetic dipole source, and 
q = u/c. We neglect the transverse polarizability of the wires assuming that they are 
thin enough. Therefore, wire polarization per unit length is solely described by axial 
polarizability a, determined from 



J 

and given by [261 1131 [H] 

-i7rql4'^ (q^R) + 



1 

a 



(4) 



(5) 



wire 

where = q^ — kl- Essential feature of Eq. ^ is the spatial dispersion of the wire 
dielectric responce, i.e., the dependence of the polarizability on the wavevector k^. 




Figure 1. Schematic illustration of the wire metamaterial with embedded light source. 
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Purcell factor for electric and magnetic dipole emission may be found via the 
imaginary part of the field, induced in the structure by the source [271 128| [25] 

/^^) =^M^^K)-P], (6) 

/(-) = ^Im[//(ro).m]. 

The field is given by a sum over Bloch waves with wavevectors = + /c^i:, similar 
as can be done for a 3D dipole lattice f2T]. The answer reads 

r_. . 3 [ Sd'k[Gu{p,)-n][Gu{-p,)-n] 
■' 2g3 7 (27r)3 - - iO " 

Here 5 = is the unit cell area, the integration over and ky is performed over the 
Brillouin zone from — vr/a to it /a, the integration range over wavevector k^ extends from 
— oo to oo and n is the unit vector oriented along the dipole. The quantity C in ([7]) is 
the interaction constant of the wires [26], defined as 

C{k)=i7rqlJ2^"'^'^Hi,'\m), (8) 

and 

Gfc(p) = ^ Go(p - P„ A:,)e'*'(''-''^) (9) 

j 

is the periodic Green function for waves with Bloch vector kj^. For electric and magnetic 
dipole problems one should substitute in ([9]) for Go electric and magnetic fields of single 
polarized wire G^^^ and Glf^\ respectively. Corresponding fields are given by 

^(p, k,) = inqlzH^'\q^p) + nk,q^pHi'\q^p) , (10) 

for electric dipole and 

G^^\p,k^) = nq^qcpHi'\q^p), (11) 

for magnetic one. Eqs. (ffl)-(fTTl) present a general explicit answer for Purcell factors. In 
practical calculations, however, the series (IH]),(l9]) should be evaluated not directly but 
using either Ewald summation [29] or Floquet-type summation [26] . The latter approach 
turns out to be more numerically efficient and yields the following result [26l [I3] 

C(fc) = -2giln^ 
a 

_ / 27iq1 sin q^^^rna gi(l - 6mfl) \ 

_ \qx,mO' cos q^^mO' — cos k^a \m\ J 



for the interaction constant, where q^^m = y/g^ — k^ — ky ^^ and ky^m = ky + 27rm/a. 
Similar technique can be used to calculate the Green function, result reads 

^k\p) ~ ^ ] [~?_L^m-^ ~ '^k^qx^mPm^ + kymkz<Smy] ■ (^3) 



Purcell effect in wire metamaterials 



5 



where 



^ ^ ^ifc^.^^ e'^"" sin q^^^x - sin q^^^{x - a) ^^^^ 

Qx 771^ COS kr^Cl cos 



_ 27r ifc„,^^ e'^"" cosg^^^a: - cosg^,^(x - a) 

e ' 

qx^ma cos fc^a - cos qx,ma 

The Green function G^™'' may be obtained from Eqs. ( fT3l) .( fT4l) using the following 
expression 

G^^^ = \xd^ + ydy + ihz) X . (15) 

Below we present numerical results for Purcell factor along with analytical answers 
in certain limiting cases along with. We first determine analytical expressions for the 
eigenmodes of the structure, Sec. [HI and then analyze Purcell factor dependence on the 
source position, Sec. 14. and on the dielectric constant of the wires. Sec. 14.21 



3. Dispersion analysis 



The spontaneous decay is due to the emission of photons, with energies found from the 
zeros of the resonant denominator in Eq. ([7]): — C = 0. Using Eqs. fj5|).f lT2|) we can 
present this denominator for ga <^ 1 and fca <^ 1 as 



-C 



where 



a 



fl2 
Hp 



fc2 



9p 



7r(l 



g2 

Hp 



27r 



+ 1 



In 



X 



71' 

+ 6. 



(16) 



(17) 



Here the wavevector x oc 1/ V^wirc — 1 characterizes the nonperfect character of the 
wires and qp is the effective plasma wavevector. The physical sense of these parameters 
is better understood when Eq. (1161) is rewritten as 



--C 
a 



4n {kl-kl){kl-kl) 



a^ql 



i.2 _|_ ^2 _ t.2 



where 



k'^ 



± 



g2 + x^ + k? 



TM 



4 



x^ 



^tm) 



g2x2 



(18) 

(19) 
(20) 



are the wavevectors of the eigenmodes of the wire medium and 



g2 - A:i - g2 . (21) 

Eigenmodes ( 1T9]) can be also obtained if the wire medium is treated as homogeneous 
medium with spatial dispersion p3], where effective dielectric constant reads 
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Figure 2. Isofrequency curves for quasi-TEM modes calculated for the different values 
of the wire dielectric constant Ewirc- Thick curves are numerical results, thin curves 
are plotted according to Eq. ([23|) . Calculation was performed at qa = O.OStt and 
R/a = 0.05. 



For large values of wire dielectric constant (£:„irc ^ 1), when x is small, Eqs. f|T9|) may 
be approximately rewritten as 



q^-ql-kl + 



(23) 



For perfect wires with e^\^c Eqs. ( I23l) reduce to 

kl = q, (TEM mode), (24) 

^2 = ^TM, (TM mode). 

The dispersionless TEM modes are specific for wire medium [13]. For q < qp they 
are the only propagating modes in the structure and thus solely control spontaneous 
emission rate. Due to the high density of TEM photonic states one can also expect large 
Purcell factor [23] . For non-perfect wires TEM and TM modes mix, as can be seen from 
Eqs. (CnD. 

Fig. [2] presents isofrequency curves of quasi-TEM modes k[^^ for different values of 
wire dielectric constant. For finite values of e^ire the TEM modes acquire hyperbolic- like 
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Figure 3. (Color online) Effective plasma frequency qp (a) and quasi- TEM mode 
wavevector (b) as functions of the wire dielectric constant. Thick black and thin red 
lines correspond to numerical calculation and Eqs. (j23p . respectively. Calculation was 
performed for Ewirc = — kwircl, qo, = O.OStt, R/a — 0.05, and xq = yo ^ a/2. 



dispersion and the absolute values of the wavevector increase according to Eq. ( |23l) . 
The growth of the wavevector for non-perfect wires is illustrated on Fig. [5^. Another 
effect of the finite dielectric constant is the decrease of the effective plasma wavevector, 
defined as the cutoff of TM waves 

ql = q'-kl{k^ = 0)^ql-x\ (25) 

Eq. (l25l) indicates, that the plasma frequency cqp decreases for non-perfect wires. 
Corresponding dependence is shown in Fig. |5)d. For small enough values of e^ire ~ {o,/ 
the plasma frequency turns to zero so the wire medium loses its properties. 

4. Purcell factor 

Now we proceed to the dependence of the Purcell factor on the dipole position within 
the unit cell of the structure and on the wire dielectric constant. 
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4--1- Effect of dipole position 

In this section we focus on perfect wires with £:„ire — oo and analyze, how the Purcell 
factor depends on the coordinate pQ = (xo,?/o) of the light source. We will first obtain 
analytical results and then compare them with numerical calculations. For perfect wires 
the integral over kz in (JTj) is determined by the residue at the wavevector corresponding 
to the TEM modes: 

m/^ |G(fc)|- JGi^a^l^ 

kl ql 

In Eq. f l26|) we have used the representation f fTSj) of the resonant denominator 1/a — C. 
For dipole close enough to the wires the dipole the Bloch Green function (j9]) is mainly 
determined by the field of the nearest wire 

Gfc(p) ^ Go(p) . (27) 

For electric and magnetic dipole this leads to 

G(^)(A;, = g) ^ -2i^e'''^ G^'\k, = q) ^ -2v^e''i' . (28) 
P P 

Eq. (125]) indicates, that for perfect wires only in-plane electric and magnetic dipoles 

couple with TEM waves. Thus, only for this orientation Purcell factor is enhanced due 

to the large density of TEM waves. Substituting Eq. fl28|) and Eq. fl26l) into Eq. ([7]) we 

obtain 

ax 

.(e) ^ ^'^X jPo ■ ^) f fe'dfc 

87r2g2 p2 J j^q2- 



Here A'max is the in-plane wavevector cutoff, steaming from the finite size of the Brillouin 
zone. We use the value -ft'max = 2-\/7r/a, so that the size of the Brillouin zone in 
cylindrical approximation stays the same, ttK^^^ = (27r/a)^. Performing the integration 
in Eq. fl29|) we obtain 

- ^''■f''' (30) 



X 



(31) 



For magnetic dipole one has 

/M(n) = /(^)(nxi), (32) 

i.e. Purcell factor for electric dipole, parallel to radius-vector Pq, is the same as for the 
in-plane magnetic dipole, perpendicular to pg. We note, that Eq. ( l32|) is the general 
relation between the Purcell factors of electric and magnetic dipole for perfect wires, 
holding at any distance from the wires. Eqs. (!30|) . (!32|) present the analytical answer for 
in-plane electric and magnetic dipoles positioned close to the wires. For vertical dipoles 
one has /i*^"* = and /i'"'* = 1 at e^ire ~^ oo, independent of the dipole coordinates. The 
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Figure 4. (Color online) Purcell factor for electric dipole as function of dipole 
position within the unit cell. Thick solid and dashed curves are numerical results 
for longitudinal and transverse dipole orientation (see geometry scheme in the inset). 
Thin lines correspond to analytical results Eq. ([501) and Eq. (|^. Calculation was 
performed for ewirc — ^ oo, qa = O.OStt and R/a — 0.05. 



general structure of Eq. fl30l) is the same as for the 3D arrays of resonant dipoles |21j . 
It consists of two factors, the first factors describes the local field enhancement. Second 
factor describes collective effect: density of states enhancement due to the TEM modes. 
It can be estimated by the order of magnitude as A^/a^, where A = 27r/g is the light 
wavelength. 

Comparison between analytical results and numerical calculation is presented in 
Fig. m Figure demonstrates high sensitivity of the Purcell factor to the position and 
orientation of the dipole. Purcell factor for electric dipole, oriented along the radius 
vector, greatly increases when the dipole approaches to the center of the wire. This 
growth is well described by Eq. ( I30l) , see thin curve in Fig. HI In case of the transverse 
orientation the Purcell factor is suppressed for small distances. Thin horizontal line in 
Fig. Hlpresents analytical answer ( l38l) for Purcell factor of the electric dipole in the lattice 
center, which is obtained in the following section and perfectly agrees with numerical 
calculation. 
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4-2. Effect of wire dielectric constant 

Fig. m shows that large Purcell factor may be obtained even for the dipole in the center 
of the unit cell. Here we focus on this case, xq = yo = a/2 and analyze the dependence of 
the Purcell factor on the dielectric constant of the wires. In this case to obtain the Green 
function it is sufficient to keep only the terms with b = in the spectral representation 
of Eq. ([9]) , which reads 



Relatively simple analytical expression for the Purcell factor can be obtained if only the 
residue at the wavevector of quasi-TEM mode with kz = ki^z is taken into account in 
Eq. ([7]). Using Eq. f llSp we perform integration over as follows: 

J 2tt (l/a) -C-iO~' An ^ ^'^^ 2k^iki - kj) ' ^ ' 

As a result, the Purcell factor is reduced to the integrals over dimensionless in-plane 

wavevector x = k±/qp: 



3g2 r ^/Y+ 



/..../(..r^Re J d.^^. (36) 





(ifmax/gp)^ 

3x'* ^ f xdx 



&^ = — Re / 



fim) _ r{m) _ Rp / 

Jy g^2^^ J 



[i+xy/^vm^ 

u 

3gp ^ f dx 



VTTxy/TT 



ux 



Here the coefficient 

ql a'^ , , 

accounts for the finite dielectric constant of the wires. The answers after the integration 
read 



.(e) 



-2- + 2 + 2v^ln ^ ^ 

a + 1 



(38) 



,,,^1,,W!_3^,1, (39, 

3o^ 2 
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Figure 5. (Color online) Purcell factor for electric and magnetic dipoles as function 
of the wire dielectric constant. Thick black, red and blue lines correspond to in- 
plane electric, vertical electric and in-plane magnetic dipole, respectively. Solid and 
dashed lines are calculated for £wiro = i£wiro and £wiro = — |£wirc|- Thin lines present 
analytical results (|38|) -(j40 |) . Calculation was performed at qa = O.OStt, R/a = 0.05, 
and xq — yo — a/2. 



In the limit of perfect wires e^ire they reduce to 

f{e) ^ f(m) _ 3 A ^ 

Jx,y J.,y 8^2^^!^^+ q2 J ^ ^^^^ 

=0, (42) 
/i™) = 1 . (43) 

Since qp ~ 1/a, the Purcell factor can be estimated as (A/a)^. There is no local field 
enhancement effect for the dipole located in the cell center, cf. Eq. ( 130|) and Eq. ( HTl) . 

Calculated dependence of electric and magnetic Purcell factor on the wire dielectric 
constant is presented in Fig. |3l The figure demonstrates, that the Purcell factor for in- 
plane electric dipole unexpectedly grows for finite dielectric constant (black curves). 
Moreover, spontaneous emission for the vertical electric dipole becomes possible and 
its rate strongly increases for smaller values of ewirc- This effect is realized both for 
superconducting wires (Ree^ire < 0, Ime^ire = 0) and for conducting wires (Re^wire = 0, 
Ime^ire > 0), cf. solid and dashed curves in Fig.El For relatively large values of l^wirel the 




Figure 6. Purcell factor for wire metamaterial, corresponding to Ref. |26j and 
operating in microwave frequency range, (a) Frequency dependence for a = 1 cm, 
R = 0.1a. (b) Dependence on lattice constant a for i? = 0.1a and v = 0.1 GHz. 



curves are well described by Eqs. fl38|) .f l39|) . The origin of this Purcell factor enhancement 
is explained by a competition of two effects: (i) dependence of the wire electric field G^^^ 
on the quasi-TEM mode wavevector fc^ and (ii) density of states dependence on kz. As 
is demonstrated by Fig. |2] and Fig. |5K, the values of for quasi-TEM modes become 
larger for non-perfect wires. Eq. ffTOj) shows, that both in-plane and axial components of 
the Green function increase with kz- The density of states for given value k± effectively 
decreases with kz, which is described by the 1/kz factor in Eq. ( I35l) . Since the Purcell 
factor is proportional to the square of the Green function times density of states, is 
still grows for smaller value of ^wire- The optimum wire of Purcell factor is reached 
at the dielectric constant |£:„ire| ~ a^/R"^, corresponding to the condition of vanishing 
effective plasma frequency gp(£:„ire), cf. Fig. [5] and Fig. [3l The situation is different 
for magnetic dipole. Contrary to the electric dipole case, the in-plane components of 
the corresponding Green function lack the pre-factor kz, see Eq. (fTTl) . Consequently, 
the Purcell factor for in-plane magnetic dipole is quenched for smaller values of wire 
dielectric constant, in agreement with Fig. |3]|^blue curves). 



5. Application to particular structures 

In this section we apply our general theory to realistic experimental structures, operating 
in microwave [2] and optical |30] frequency ranges. The structure from Ref. [2], used for 
subwalength transmission of images, is characterized with the period a = 1 cm and wire 
radius R = 1 mm. Asymptotic answer Eq. (HTl) for perfectly conducting wires is valid 
for microwaves. This answer may be explicitly written as 
r = f- = 37rc^ ln[l + 7r/3 + 21n(a/27rfi)] 

J^,y J-,y 4^2^2 ln(a/27ri?) + 7r/6 ' ^ ^ 

First factor in Eq. fj44|) mainly determines the dependence of the Purcell factor on 
structure parameters, while the second factor weakly depends on the relative thickness of 
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Figure 7. Purcell factor for nanowire metamaterial, corresponding to Ref. [30]. (a) 
Energy dependence of the dielectric constant of the wires, (b) Energy dependence of 
the Purcell factor for horizontal (black solid curve) and vertical (red dashed curve) 
electric dipole. Calculated for a = 60 nm, R — a/4, xo = ^ a/2. Decay rate is 
normalized to the value in bulk matrix. 



the wires R/a. Dependence of the Purcell factor on the the source frequency and on the 
lattice constant is shown in Fig. EJ^a) and Fig. EJ^b), respectively. Purcell factor strongly 
increases for smaller frequencies and smaller lattice periods. Calculation demonstrates 
that the wire medium allows to achieve high values of Purcell factor in the whole 
microwave spectral range. 

In the optical frequency range wire metamaterial may be realized as an array of gold 
nanowires in alumina [301 E]. Purcell factor calculated for this structure is presented 
in Fig. [3 Reahstic energy dependence of the gold dielectric constant of nanowires 
from Ref. [30] has been taken into account. We have used the value ^out = 2.56 for 
the dielectric constant of the alumina matrix. To account for the dielectric constant 
of the matrix, different from unity, wire dielectric constant in Eq. (j2]) was scaled as 
f^wire — ^ ^wire/f^out, which yiclds the spontaneous decay enhancement with respect to the 
bulk matrix. Fig. [7]shows, that for axial dipole orientation a flat maximum in the Purcell 
factor with / ~ 20 is reached at the energy £^ ~ 1.8 eV. This roughly corresponds to the 
optimal dielectric constant of the wires, revealed in Fig. |3l For in-plane dipole the Purcell 
factor monotonously decays with photon energy due to decreasing photonic density of 
states ~ l/(ga)^, although its value stays above 10 for all energies below 2 eV. Several 
issues should be noted regarding the relevance of Fig. [7] to real experimental structure. 
First, our theory has been derived for thin wires, neglecting their transverse polarization. 
This assumption certainly fails at large energies. Second, the real nanorods have the 
lengths in order of hundreds of nm, while in theory they are assumed infinite. Third, 
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the Purcell factor calculated according to Eq. (E]) corresponds to the total decay rate 
caused by electromagnetic interaction with the medium. It is determined by a sum of 
the rate of photon radiation in the far field and the rate, with which photons are radiated 
and then reabsorbed due to the dielectric losses [25]. Thus, Eq. ([6]) overestimates the 
enhancement of radiation efficiency. Consequently, curves in Fig. [7] should be considered 
as an upper boundary of the Purcell factor in gold nanorod metamaterial, rather than 
as a rigorous modeling. Still, they indicate that relatively high Purcell factor may be 
achieved in the visible spectral range. 

6. Conclusions 

We have developed a general analytical theory of spontaneous emission of both electric 
and magnetic dipole sources in wire metamaterials. Our theory goes beyond the effective 
medium approximation and fully accounts for the discreteness of the structure. We have 
analyzed the dependence of the Purcell factor on the dipole position within the lattice 
unit cell as well as on wire dielectric constant. We have demonstrated that the Purcell 
factor can be greatly enhanced due to the large density of states of TEM modes, and 
its values are of the order of a square of the ratio of the light wavelength and lattice 
constant. The Purcell factor is also very sensitive to the position and orientation of 
the dipole source, and it may increase due to the local field effect when the dipole 
approaches the wires. Counterintuitively, the spontaneous emission rate for the electric 
dipole emission grows when the wires are not perfectly conducting. We have found 
optimal value of the wire dielectric constant that maximizes the Purcell factor. We 
have demonstrated the possibility of broadband spontaneous decay rate enhancement 
in realistic wire metamaterials operating in both microwave and optical spectral range. 
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